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Abstract 

The RHIC program was intended to identify and study the quark-gluon plasma formed in the 
collision of heavy nuclei. The discovery of the "perfect liquid" is an essential step towards the 
understanding of the medium formed in these collisions. Much of data relevant to this was 
provided by the study of "soft" observables, which involve many particles of low momentum 
produced in nearly every event, rather than high momentum particles produced in rare events. 
The main results related to soft physics at RHIC are discussed, as well as their implications for 
the physics of the LHC heavy ion program. 
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1. RHIC physics in a nutshell: The "perfect liquid" 

The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory turned 
on in 2000 and has taken data over a wide range of energies (19.6 to 200 GeV per nucleon- 
nucleon collision) and collision systems (protons to gold). The expectation from lattice 
QCD calculations (an example shown in the left panel of Fig. [I] from Ref. pQ) was that 
colliding nuclei at relativistic energies would form a hot, dense system where the degrees 
of freedom would no longer be the hadrons measured in the final state, but rather their 
quark and gluon constituents. However, the data collected by the four experiments at 
RHIC - two large, and two small, all covering both unique and overlapping ranges in 
phase space - arrived at the surprising conclusion that the system formed in the colli- 
sions was a "perfect liquid" [2J. This is a non-trivial observation, given that asymptotic 
freedom implied that interactions between the quarks and gluons should become weaker 
at higher energies [3J. Instead, the system appears to flow collectively as if the interac- 
tions between the relevant degrees of freedom are exceedingly strong. The latter property 
can be quantified by the viscosity, which is large for an non-interacting gas, and zero at 
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Fig. 1. (left) Lattice QCD calculations showing the transition near T c = 170 MeV from Ref. |T . (right) 
Viscosity scaled by entropy density as a function of temperature for three common liquids (He, N, H2O), 
normalized such that the viscosity bound is at unity. 

infinite coupling. Early estimates of the viscosity of the medium formed at RHIC (dis- 
cussed below) suggest that it is within a factor of 2 from the so-called "viscosity bound" 
predicted using the AdS/CFT correspondence [4], shown compared to normal fluids in 
the right panel of Fig. [T| This may well be the first prediction from string theory to 
be validated by experiment, a major development for both heavy ion physics and string 
theory. 

This work will discuss the most important results of "soft" physics at RHIC, which 
involve the characteristics of multiparticle production at low momentum, rather than the 
rare production of high momentum particles. Soft physics is of great practical value, as it 
is available to experiments right after a machine starts. And while it is not yet amenable 
to perturbative approaches, it provides a handle on a sector of QCD of great interest for 
understanding the medium produced in RHIC collisions. 

2. Evidence for the perfect liquid from soft observables 

The discovery of the perfect fluid at RHIC was established by a series of inferences 
from the final state back in time to the initial state. The spectrum of final state pions 
is manifestly blackbody in shape, as one can see in distributions of identified particles 
(e.g. Ref. |S]). Integrating the blackbody spectrum for the different hadronic mass states 
(with proper Bose and Fermi distributions), and taking ratios to factor out the emission 
volume, gives "thermal model" [B] predictions for ratios of total yields of different hadron 
states. Comparison to measurements, an example of which is shown in Fig. [2] from the 
STAR experiment [7] , shows an excellent agreement and requires only three parameters: 
a chemical frcczcout temperature (T c h), a baryochemical potential (/is) reflecting the 
excess of baryons in the initial state, and a strangeness suppression factor (which is 
found to be nearly unity at RHIC). These data suggest that the system was thermalized 
at the time of freezeout and was thus at least T = 170 MeV during its evolution. 

The success of describing the data from such a wide variety of hadron states validates 
the approach begun by Rolf Hagedorn in the 1960's, who observed the thermal slopes in 
hadron momentum spectra and postulated that they were emitted from a thermalized 
state [8]. The observation of the same temperature in all systems led him to also propose 
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Fig. 2. (left) STAR data on hadron ratios fit to a thermal model, from Ref. 7 . (right) The elliptic flow 
parameter V2 vs. N par t for y/sjvjv = 200 GeV from PHOBOS compared with a hydrodynamical model. 

a "limiting temperature" arising from an exponentially rising mass spectrum (dN/dM oc 
exp(Af/Ty)). The rising mass spectrum was a prediction borne out by forty years of 
subsequent measurements, making surpassing the Hagedorn temperature Tjj ~ 170 MeV 
that much more interesting for heavy ion experiments forming hot and dense matter. 
Lattice QCD calculations appear to have no limiting temperature, reflecting the fact 
that they find a discontinuity in the number of degrees of freedom, e.g. as quantified by 
e/T 4 shown as a function of T in the right panel of Fig[TJ at a temperature quite similar 
toT H p. 

Beyond being thermal, heavy ion collisions show a collective flow characteristic of a 
liquid. At all energies and system sizes, it has been observed that the event-by-event 
angular distributions are not isotropic in azimuthal angle. An "event plane" can be 
estimated from the produced particles, defined as angle *S?r of the short principal axis 
of the particle angles. The azimuthal distribution relative to this event plane is found to 
show a strong cos(2[0— ^r]) dependence. This is especially pronounced in more peripheral 
(i.e. lower multiplicity) collisions, where the overlap of the nuclei is shaped like an almond, 
relative to central (i.e. higher multiplicity) collisions where the overlap is essentially 
isotropic. This leads to a characterization of the event-by-event angular distributions in 

terms of its Fourier coefficients [9] = 1 + 2«i cos(0 - Wr) + 2v 2 cos(2[0 - V R ]) H . 

PHOBOS data on V2 is shown in Fig. [2] and compared to a hydrodynamical calculation 
tuned to the most central data and extrapolated to larger impact parameters using a 
Glauber model [TO] . 

Comparing the relevant energy and space-time scales implied by the success of the 
hydrodynamical models, the matter at RHIC is formed under quite extreme conditions. 
The formation time needed for the the hydrodynamic calculations is tq — 0.6 fm/c, or 
approximately 2 yoctoseconds (10 -24 ) [TT]. This number is far smaller than the time 
taken a massless particle to traverse the radius of a hadron (r ~ 1 fm/c) [12]. The same 
calculations determine that the energy density needed to match the data is around e ~ 30 
GeV/fm 3 , about 60 times the density of a nucleon in its rest frame, ejy ~ 500 MeV/fm 3 . 
It should be noted that these estimates do not preclude even higher energy densities at 
even earlier times. 

The agreement of hydrodynamic models with RHIC collisions is an important, but 
surprising, development. Hydrodynamics is not merely "another" model, but describes a 
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medium with particular properties, difficult to generate by the rescattering of partons or 
hadrons. It is a locally equilibrated system, with a well-defined temperature at any given 
space-time location. The agreement with data, particularly the detailed flow of heavy 
particles, seems to require a "lattice-like" equation of state, with a speed of sound that 
drops rapidly at the hadronization temperature. With an energy density at least 60 times 
the energy density of a proton, the system resists a description as interacting hadrons, 
as their wave functions would all be overlapping and thus indistinguishable as individual 
particle states. 

This same energy density implies a temperature about 2-3 times the critical tempera- 
ture found on the lattice, well above the Hagedorn temperature. Exceeding the Hagedorn 
temperature is a key piece of evidence required for any system claiming to have degrees 
of freedom comprised of quarks and gluons. The recent PHENIX data, shown in the left 
panel of Fig. '3, show a spectrum of direct virtual photons from heavy ion collisions, com- 
pared with a similar sample from proton-proton collisions |13j . There is a notable excess 
in the lowest px bins which can be fit by a simple exponential of slope T ~ 220 MeV. 
These data do not fully establish the temperature of the system at the very early times, 
as the measured spectrum is a convolution of spectra emitted during the whole collision 
evolution. Still, they directly establish that the emitted photons were hotter than the 
Hagedorn temperature at some point in the collision evolution. This suggests that the 
medium may well be the QGP, although this result per se does not point directly to quark 
and gluon degrees of freedom. Interestingly, comparisons of data for pions and protons 
shown in the right panel of Fig. [3] require a sharper transition than actually observed in 
lattice calculations [3"fj] . 



3. The edge of liquidity: Can we turn off the perfect liquid? 

The detailed characterization of the medium formed at RHIC requires a systematic 
variation of the initial conditions to measure its effects on the final state particles. Such 
scans have been performed in other systems, e.g. ultra-cold atomic gases [14] . Since soft 
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Fig. 4. (left) Extended longitudinal scaling observed for inclusive charged particle distributions in 
Au+Au collisions over a wide range of energies 1191 , (right) Extended longitudinal scaling seen in e+e - 
reactions compared to heavy ion data (figure from Ref. 1201 ), 



physics involves the bulk production of particles from a thermalized system undergoing 
collective flow, one would want to vary the temperature, density and material properties 
in such a way to quantitatively extract the equation of state. Unfortunately, the only 
tools available for heavy ion physics are colliding ion beams, for which one can only 
control the beam energy and nuclear size. It also possible to vary the centrality of the 
collision, which affects the shape of the overlap region. 

The most fundamental question addressed in all of these experimental studies is if the 
system is sufficiently thermalized in the initial state, such that it can properly be called 
a state of matter. Once that is established, a whole series of questions arises, each of 
which addresses fundamental problems in strong interaction physics. 

- Does the system thermalize everywhere, or just in limited regions of phase space? 

- What are the conditions (size, energy) for thcrmalization? 

- How rapidly is thcrmalization achieved? 

- What degrees of freedom reach thermalization (partons, hadrons, or something else)? 
If it is assumed that the top-energy, most-central RHIC collisions are in fact thermalized 
(as suggested by their agreement with hydrodynamical models), it should be possible to 
vary the experimental parameters and "turn off" the creation of quark gluon plasma. 

Lowering the beam energy would seem to be a natural way to create a system that does 
not form a quark gluon plasma. However, it is essential to check if this leads to particle 
production different in some essential way than at higher energies. The inclusive charged 
particle multiplicity is a useful measurement for this, as it is a "global" observable, which 
integrates over time, space, and various degrees of freedom. The number of particles 
emitted essentially counts the degrees of freedom available to the system, which should 
be linear with the produced entropy in a fully thermalized system |16j . One empirically 
observed feature that should shed light on this is "extended longitudinal scaling" 17 \. It 
has been observed in proton-proton collisions that the pseudorapidity density {dN c h / drj) 
of inclusive charged particle production is energy-independent when viewed in a frame 
where one or other of the incoming particles is at rest |18j . This is done by using the 
kinematic variable rj' = r\—ybearm where ybeam is the rapidity of one of the beams. Results 
for dN/drf are shown for Au+Au collisions at four RHIC energies [T^ in the left panel 
of Fig. [4] where longitudinal scaling is clearly observed. The persistence of this scaling 
over a factor of ten in energy suggests that no major changes in the particle production 



5 



(a) 

4 



: | n 200GeV Au-Au /r ^) 

15 1 T ■ 62.4GeV Au-Au If 

— t-OOftmnc ^ 2D0GeV Cu-Cu ** 

E - ] • 62.4GeV Cu-Cu If 



(c) 



200 
N parl 



occurs over this range. It also suggests that physics at rj = is not obviously from a 
different origin than physics in the forward direction. 

What is surprising is that the same phenom- 
ena is common to all systems, from heavy ions 
to proton-proton, and even to e + e~ annihilation. 

The latter system is shown in the right panel of >"„.„„ i l|f | J * 

Fig. [4] which shows that the phenomenon is also f 
not dependent on the overall system size [15] . It 
also shows an overlay of three Au+Au energies [20 , 
showing that the overall magnitude of multiplic- 
ity in heavy ion collisions and e + e~ reactions is 
essentially the same, provided the Au+Au data 

is scaled down by the number of participant pairs o. 5 - f . 

(Np ar t/2). While this similarity between heavy ions : §t*n . ■ * * 

and e + e _ has been reported for several years |21j . o.*~ 7 

it has not been fully explained. ^ a . 3 - 

PHOBOS data in Ref. [H] show that the el- «\,L 

liptic flow parameter V2 also obeys extended Ion- * „ , L ^ * 

gitudinal scaling. This is surprising if the initial f 

. ... 100 

conditions at midrapidity are completely differ- 
ent to that at forward rapidities, but not if the 

initial entropy density and geometry are the pri- ^OBOSdJ&m 4 ' ^ 
mary determinants of elliptic flow, as should be 

the case if the viscosity is small. It has been shown by various experiments that the 
elliptic flow, normalized by the initial eccentricity, scales with the transverse particle 
density dN/dy/S, where S is the area overlap of the initial collision [24J. However, this 
scaling breaks down when comparing Au+Au and Cu+Cu reactions if the "standard" 
eccentricity is used: e st d — (<7y — cr x)/(°i' + CT jf)> wnere tT x( cr y) is the variance in the 
direction along (perpendicular to) the reaction plane. Figj5]shows that scaling is restored 
by using an eccentricity defined by the distribution of the participants themselves [23 : 
e P art = (\/ ( a Y ~ a x) 2 + 4(°xy ) 2 )/(°x + °y)- This suggests that the geometrical con- 
figuration of the participants is "frozen in" immediately, consistent with the previous 
estimates of To or perhaps even shorter times. 

The predictions from AdS/CFT for shear viscosity [3] have made measurements of 
observables sensitive to viscous effects a high priority. Currently, there are two approaches 
for making quantitative estimates of r\j s based on experimental data. One is based on 
direct measurement of non-equilibrium processes, e.g. the behavior of charm quarks which 
are produced in the earliest stages. Model calculations arc used to relate the observed 
magnitude of V2 with simultaneous measurements of suppression as a function of pt- 
PHENIX results [37], shown in Fig(6] find good agreement between the data and models 
assuming that 77/s ~ (1 — 2)/47r. The other method is to systematically include viscous 
corrections into the hydrodynamic evolution, a technically challenging task which is being 
addressed by several competing teams [2S] . An example of this from Ref. [25] is shown 
in the right panel of Fig. [6] and compared to PHOBOS data. The main uncertainty in 
these calculations at present is the initial state (e.g. whether the energy is distributed 
according to Glauber models or according to the Color Glass Condensate) . However, the 
effect on rj/s is only about a factor of two, which suggests that RHIC is quite close to 
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Fig. 6. (left) PHENIX data Ref. 1271 on charm suppression and elliptic flow for Au+Au collisions, (right) 
PHOBOS data compared to a set of calculations using the viscous hydro code from Ref. 1201 . 
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Fig. 7. (left) PHENIX data on V2/n q (from Ref. I28I . divided by a reference) vs. KEj-/n q for various 
hadrons (right) the same quantity for a hydrodynamic calculation, both from Ref. |29| . 



the viscosity bound. 



4. Degrees of freedom: What is the perfect liquid? 

If the near-perfect fluid description is relevant over most of the evolution, it is essential 
to determine the fluid constituents, and how they become a thermalized collective state of 
matter. Nothing in the data discussed so far uniquely identifies which degrees of freedom 
are able to achieve this. The most natural assumption would be that the early stages 
are dominated by the dynamics of free quarks and gluons, or at least the dynamics of 
quark and gluon fields that are studied using Lattice QCD. However, as shown above 
attempts to model the existing data on vs. pt for identified particles find that a strong 
first-order phase transition needs to be put in by hand [30 . Existing lattice calculations, 
shown in Fig.[3]do not provide sufficient "softening" of the equation of state to reproduce 
the heavier particles which are most sensitive to the speed of sound. 

Another way to look for quark degrees of freedom is via constituent quark (n q ) scaling. 
The i>2 data have been studied for many different particles species, which have different 
mass and quark content. PHENIX data [28 show that all of the available data on v-i 
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Fig. 8. (left) Data on mid-rapidity multiplicity in p+p and A+A, with predictions for the LHC energies, 
from Ref. 1351 . (right) Coverage of LHC experiments in the forward direction, presented in the rest frame 
of one of the projectiles. 

vs. pt lie near one another when plotted as a function of V2/n q on the Y axis, where n q 
is the number of valence quarks and anti-quarks in the hadron, and KEx/riq on the X 
axis, where KEj> — — m for each hadron of mass m. This suggests a scenario where 
freezeout occurs by the recombination of constituent quarks, assumed to have a mass 
of ~ m/n q and the appropriate quantum numbers for each hadron. However, Ref. |29l 
argues that the presence of good quasiparticles, where the width is small relative to the 
mass, should induce viscous effects which could exceed the estimates of viscosity using 
the data in Fig. [6] And yet, the V2/rig data scaled by a reference derived from the kaon 
flow (shown in the left panel Fig. PFI show that the flow data fit into the n q scaling 
pattern, while hydro calculations (shown in the right panel) unsurprisingly do not |29j . 
These data imply that it remains an open question how to harmonize the perfect fluid 
and constituent quark scenarios. If one tries to recapture the fluid behavior by allowing 
the large quasiparticle widths, it is not clear at which point the quasiparticles no longer 
act as the relevant degrees of freedom. 



5. The future: RHIC II and the LHC 



The perfect liquid is now the paradigm that will be extensively tested by the next 
generation of heavy ion experiments. These will take place at an upgraded RHIC II, with 
a factor of ten increase in luminosity, and at the upcoming LHC Pb+Pb program, with a 
factor of nearly thirty increase in center-of-mass energy. Probing the transport properties 
of the system with jets and heavy quarks will be the focus of the next generation of RHIC 
experiments [31] ■ New silicon detectors in PHENIX and STAR are being developed to 
measure charmed particles by means of displaced decay vertices. The LHC will explore 
both hard and soft physics in a new energy regime, with the three large LHC experiments 
- ALICE, ATLAS and CMS - all with active heavy ion programs |32|33|34| . Hard probes 
will be especially powerful, due to the increased rate of jets, photons, and quarkonia 
expected from perturbative QCD calculations. 

However, soft physics will also lead to great strides in our understanding of the strongly- 
coupled medium. Furthermore, the results will come out quickly, and be able to imme- 
diately address the relevance of models and extrapolations of lower energy data. A key 
observable is the charged particle multiplicity, for which data and three predictions are 
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shown in the left panel of Fig.[8j from Ref. [35] • The coverage of ATLAS in the rf variable, 
relevant for testing extended longitudinal scaling, is shown in the right panel. Whether or 
not the data agree with current expectations, there will surely be new insights on the na- 
ture of the perfect fluid. Of particular interest is the possibility that the system continues 
to thermalize very early, something which will not occur if perturbative physics becomes 
more important at higher energies. The diversity of experiments and technologies will 
guarantee that the available phase space is fully covered for a variety of observables, 
making the early days of LHC physics in both p+p and Pb+Pb particularly exciting. 
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